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ABSTRACT 

The Vela pulsar represents a distinct group of 7-ray pulsars. Fermi 7-ray observations reveal that it has 
two sharp peaks (PI and P2) in the light curve with a phase separation of 0.42 and a third peak (P3) in the 
bridge. The location and intensity of P3 are energy-dependent. We use the 3D magnetospheric model for the 
annular gap and core gap to simulate the 7-ray light curves, phase-averaged and phase-resolved spectra. We 
found that the acceleration electric field along a field line in the annular gap region decreases with heights. The 
emission at high energy GeV band is originated from the curvature radiation of accelerated primary particles, 
while the synchrotron radiation from secondary particles have some contributions to low energy 7-ray band 
(0.1 — 0.3 GeV). The 7-ray light curve peaks PI and P2 are generated in the annular gap region near the 
altitude of null charge surface, whereas P3 and the bridge emission is generated in the core gap region. The 
intensity and location of P3 at different energy bands depend on the emission altitudes. The radio emission 
from the Vela pulsar should be generated in a high-altitude narrow regions of the annular gap, which leads to a 
radio phase lag of ^ 0. 13 prior to the first 7-ray peak. 

Subject headings: Pulsars: general - Gamma rays: stars - radiation mechanisms: non-thermal - Pulsars: 
individual (PSR J0835-4510) 



L INTRODUCTION 

The Vela pulsar is the brightest point source in the 7- 
ray sky. The Vela pulsar at a distance of d = 287l]^7pc 
(iDodson et alj|2^03i) has a spin period of P = 89.3 ms, char- 
acteristic age Tc = llkyr, magnetic field B = 3.38 x 10^^ G, 
and the rotational energ y loss rate E = 6.9 x lO'^^ergs"^ 
(iManchester et al.ll2005h . It radiates multi-waveband pulsed 
emission from radio to 7-ray which enables us to get consid- 
erable insights of the magnetosphere activities. High energy 
7-ray emission ta kes away a significant f r action of the spin - 
down luminosity ( [Thompson et al.l [19991 : iThompsonI 1200 Ih . 
The pulsed 7-ray emission from the Vela pulsar was detected 
by many instruments, e.g. SAS 2 (iThompson et al.l [19751) . 
COS B (Grenier et an iI988h . the Energetic Gamma Ray Ex- 
periment Telescope (EGRET, Kanbach et al. 1994; Fierro 
et al. 1998), Astro-rivelatore Gamma a I mmagini LEggero 
(AGIL E, Pellizzoni et al. 2009) and Fenni jAbdo et alj|2009l 
l20I0bl) . The 7-ray profile has two main sharp peaks (PI and 
P2) and a third peak (P3) in the bridge. The location and in- 
tensity of P3 as well as the peak ratio (P1/P2) vary with energy 
jAbdo et alj|2010blc[) . Because of the large E, the Vela pul- 
sar has a strong wind nebulae, fro m which the unpulsed 7 -ray 
photons was detected by AGILE (iPellizzoni et alJl20IOh and 
Fenni (Abd_o_et al. 2010al). 

Theories for non-thermal high energy emission of pulsars 
are significantly constrained by sensitive 7-ray observations 
by the Fenni telescope. Four physical or geometrical mag- 
netospheric models have previously been proposed to ex- 
plain pulsed 7-ray emission of pul sars: the polar cap model 
jDaughertv & Harding 1 19941 1 1 996h in which the emission re- 
gion is generated near the n eu tron st ar surface, the out er gap 
model (Cheng et al. 198 6a|[bl: iRom ani & Yadigaroglul 
Zhang & Chena, 1997; , Cheng et alJi200(l: .Zhang et all 
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in which the emission region is generated near the light 
cylinder, the two-pole caustic model or the s l ot gap model 
(IDvks & RudaS l200l IMushmov & Harding! l2003l l200l 
[Harding et al.1 120081) in which the emission region is gener- 
ated alo ng the last open fie ld lines, and the annular gap model 
dOiao e t al. 2004a b, 2007: iDu et al]|2010l) in which the emis- 
sion is generated near the null charge surface. The distin- 
guishing features of these models are different acceleration 
region for primary partic les and possible mechanisms to radi- 
ate high energy photons. iRomani & Yadigaroglul(II995l) mod- 
eled the 7-ray and radio light curves for the Vela pulsar with 
a larger viewing angle ^ 79°). In their outer gap model, 
the two 7-ray peaks are generated from the outer gap of one 
pole, whereas the radio emission is radiated from the other 
pole. However, the correlation of high energy X -ray em ission 
and the radio pulse shown by Lommen et all (120071) is not 
consistent with this picture. iDvks & RudakI (120031) used the 
two-pole caustic model to simulate the 7-ra y light curve for 
the Vela pulsar, wh i ch wa s further revised by lYu et al.l (120091) 
and lFang & Zhang! (120 1 Ol) to explain the details of Fenni GeV 
light curves. A bump appears in the bridge in the model for 
a large inclination angle, but the width and the location of P3 
were not well modeled yet. 

In this paper, we focus on the 7-ray light curves at diffemt 
bands and spectra of the Vela pulsar. In § 2, we introduce the 
annular gap and core gap, and calculate the acceleration elec- 
tric field in the annular gap. In § 3, we model the multi-band 
light curves using the annular gap model together with a core 
gap. The radio emission region is identified and the radio lag 
prior to the first 7-ray peak is explained. To model the Vela 
pulsar spectrum, we also calculate the 7-ray phase-averaged 
and phase-resolved spectra of both synchro-curvature radia- 
tion from the primary particles and sychrotron radiation from 
the secondary particles. Conclusions and discussions are pre- 
sented in § 4. 
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2. THE ANNULAR GAP AND CORE GAP 
2.1. Formation of the Annular Gap and the Core Gap 
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Fig. 1 . — Sketch for the annular gap and the core gap for particle acceler- 
ation. CFL stands for the critical field line across the intersection of the null 
charge suiface and the light cylinder, NCS for the null charge suiface, and 
LOF for the last open field line. P is the peak emission spot at a field fine in 
the annular gap region, which is loc ated between CFL and LOF. This figure 
is taken and revised from Figure 1 of lOiao et al.ld2007l) . 

The open field line region of pulsar magnetosphere can be 
divided into two parts by the critical field lines (see Figure[T]l. 
The core region near the magnetic axis is defined by the crit- 
ical field lines. The annular region is located between the 
critical field lines and the last open field lines. For an anti- 
parallel rotator the radius of the core gap (rcoro) and the full 
polar cap region (rp) are rcorc = (2/3)'^/^i?(rii?/c)^/^ and 
rp = R{QR/cy/^, respectively (iRuderman & Sutherland! 
11975b . where R is the neutron star radius, Q is the angular 
velocity {^l = 2tt/P, P is the pulsar spin period). The radius 
of the annular polar region therefore is rann = ''p ^ ''coic = 
0.26i?(rii?/c)^/^. It is larger for pulsars with smaller spin 
periods. 

The annular acceleration region is negligible for older long 
period pulsars, but very important for pulsars with a small pe- 
riod, e.g., millisecond pulsars and young pulsars. It extends 
from the pulsar surface to the null charge surface or even be- 
yond it (see Figure[TJ. The annular gap has a sufficient thick- 
ness of trans-field lines and a wide altitude range for parti- 
cle acceleration. In the annular gap model, the high energy 
emis sion is generat ed in the vicinity of the null charge sur- 
face dP u et al. 201tf). This leads to a fan-beam 7-ray emission 
jOiao et al. 2007). The radiation from both the co re gap and 
the ann ular gap can be observed by one observer (lOiao et al.l 
l20041j) if the inclination angle and the viewing angle are suit- 
able. 

2.2. Acceleration Electric Field 

The charged particles can not co-rotate with a neutron star 
near the light cylinder and must escape from the magneto- 
sphere. If particles escape near the light cylinder, these par- 
ticles have to be generated and move out from the inner re- 
gion to the outer region. This dynamic process is always tak- 
ing place, and a huge acceleration electric field exists in the 
magnetosphere. To keep the whole system charge-free, the 
neutron star surface must supply the charged particles to the 
magnetosphere. 

The annular gap and the core gap have particles with oppo- 
site sign flowing, which can lead to the circuit closure in the 
whole magnetosphere. The potential along the closed fiel d 
lines and the critical field lines are different (IXu et al.ll2066h . 
The parallel electric fields {E0 in the annular gap and core 
gap re gions are opposite, as has been discussed by ISturrockl 
(11971b . As a result, E'n vanishes at the boundary (i.e. the crit- 
ical field lines) between the annular and the core regions and 
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Fig. 2. — For the Vela pulsar, the acceleration electric field is calculated 
along the field line marked as the dotted line with 'P' in Figure [T] with a 
magnetic azimuthal ip = 0° and a magnetic inclination angle a = 70° in 
the annular gap model. 



also along the closed field lines. The positive and the negative 
charges are accelerated from the core and the annular regions, 
respectively. 

We now consider a tiny magnetic tube in the annular gap 
region. We assume that the particles flow out at a radial 
distance about rout ~ ^lc — 4.3 10^ km, and that 
the charge density of flowing-out particles Ph{rout) is equal 
to the local Goldreich-Ju lian (GJ) charge density /9gj(rout) 
dOoldreich & JulianI 19691) at a radial distance of rout- For any 
heights r < rout, Phif) < Pgjir)- The acceleration electric 
field therefore exists along the field line, and cannot vanish 
until approaching the height of rout- 

For a static dipole magnetic field, the field components can 
be described as = ^^^^n^ and Bg = ^^n^, here 9 
is the zenith angle in magnetic field coordinate, and Bq is the 
surface magnetic field. Thus the magnetic field at a height r is 
B{r) = Bo^VM^m^ In the co-rotating frame, Poisson's 
equation is 

V-E = 4^(pb-Pgj)- (1) 

Because of the conservation laws of the particle number and 
magnetic flux in the magnetic flux tube, the difference be- 
tween the flowing charge density and local GJ charge density 
at the radius r can be written as 

ilB(r) 

Ph{r) - Pgj(r) = 2^*'*^°^''°"* ~ ^^-^ 

where Vl = 2tt/P is the angular velocity, P is the rotation 
period, and (and Cout) are the angle between the rotational 
axis and the B field direction at r (and rout)- Wang et al. 
(2006) found 

cos C = COS a cos 0^ — sin a sin 0^ cos "0, (3) 

where ip and are the azimuthal angle and the tangent angle 
(half beam angle) in the magnetic field coordinate, respec- 
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tively. Combining equations ([T]i, (|2]i and ([3]l, we obtain 



V-E = -- 



Substituting cot 



- Vocos^^ + l (cos (out — COS C) • (4) 



2 cot^ ( 



(IOiao&Linlll998h and ds 



M ~ 3 cote 

y'(7YM)2 + (ci7^ into equations ([3j and (|4|, We can solve the 
equation for V • E, and calculate the electric field E'n along a 
magnetic filed line for ip ~ 0°, as shown in Figure Ofor the 
Vela pulsar The electric field is huge in the inner region of 
annular gap and drops quickly when r ^ Rhc- 

3. MODELING THE Fermi 7-RAY PROFILES AND 
SPECTRA OF THE VELA PULSAR 

We reprocessed the Fermi data to obtain the multi-band 
light curves in the following steps: (1) Limited by the tim- 
ing solution for the Vela pulsafl from the Fermi Science Sup- 
port Center (FSSC), we reprocessed the original data ob- 
served from 2008 August 4 to 2009 July 2. (2) We selected 
photons of 0.1-300GeV in the "Diffuse" event class, within 
a radius of 2° of the Vela pulsar position (RA= 128.55°, 
DEC= -45.75 °) and the z enith angle s maller than 105°. 
(3) As done by lAbdo et al.l (120091 l2010b llB), we used "fse- 
lect" to select photons of energy Eq^v within an angle of < 
max[1.6 — 3 log]^g(iJGov), 1-31 degrees from the pulsar posi- 
tion. (4) Using the tempo2 (iHobbs et al.ll2006t lEdwards et all 
l2006h with the Fermi plug-in, we obtained the rotational phase 
for each photon. (5) Finally we obtained the multi-band 7-ray 
light curves with 256 bins, as presented in Figure |4](red solid 
lines). Two sharp peaks have a phase separation of 5(f) ^ 0.42. 
The ratio of P2/P1 increases with energy. A third broad peak 
appears in the bridge emission. The intensity and phase loca- 
tion of P3 vary with energy. 

These observed features challenge all current high energy 
emission models. A convincing model with reasonable input 
parameters for magnetic inclination angle a and viewing an- 
gle C should produce multi-band light curves of the Vela pul- 
sar and explain the energy-dependent location of P3 as well 
as the ratio of P2/P1. 

3.1. Geometric Modeling the Light Curves 

Model parameters for both the annular gap and core gap 
of the Vela pulsar should be adjusted for the particle accel- 
eration regions where the 7-ray emission are generated. The 
framework of the annular gap model as well as the coordinate 
details have been presented in lDu et all (120101) . which can be 
used for simulation of the multi-band 7-ray light curves of 
the Vela pulsar In this paper, we added the simulations for 
the core gap to explain P3 and bridge emission. We adopted 
the inclination angle of a = 70° and the viewing angle 
C = 64° wh ich were obtained from the X-ray torus fitting 
dNg & Rom ani 2008). The modeling was done as follows. 

1 . We first separate the polar cap region into the annular and 
core gap regions by the critical field line. Then, we use the 
so-called "open volume coordinates" (rove 7 i's) to label the 
open field lines for the annular gap and core gap, respectively. 
Here rove is the normalized magnetic colatitude and V-'s is 
the magnetic azimuthal. We define ij^s — for the plane of 
the magnetic axis and the spin axis, shown in Figure 1 . For 
the annular gap, we define the inner rim rove, AG = for 
the critical field lines and the outer rim rove. AG = 1 for 

' http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ 
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Fig. 3. — The height of the emission peak rp of the Vela pulsar in the 
annular gap model and the height of null charge surface rN, calculated with 
an incliantion angle a = 70°, n = 0.7 and A = 0.9. Note that rp and r^j 
are symmetric around the magnetic axis in the magnetic frame. The projected 
Tp is always within the light cylinder. We define i/ig = 0° for the median 
between the magnetic axis and the equator in the plane of the spin axis and 
magnetic axis. 



the last open field Unes; while for the core gap, we define the 
outer rim rove, eG = 1 for the critical field lines and the 
inner rim rove, ee = for the magnetic axis. We also divide 
both the annular gap (0 < rove, AG ^ 1) and the core gap 
(0.1 < rove, ee ^ 1) into 40 rings for calculation. 

2. Rather following the conventional assumption of the 
uniform emissivity along an open field line when modeling 
the light curve s (Dvks & Rudak 2003; Harding eTaI1l2008t 
iFang & Zhang||2010l) . for both the annular gap and the core 
gap, we assume that the 7-ray emissivities /(6's, "08 ) along one 
open field line have a Gaussian distribution, i.e.. 



2al 



(5) 

here 6s is the magnetic colatitude of a spot on a field line, 
V's is the magnetic azimuthal of this field line, C{9s,ips) = 

JH" -^/r^ + (dr/dd)^ dO is the arc length of the emission 
point on each field line counted from the pulsar center, cta 
is a length scale for the emission region on each open field 
line in the annular gap or the core gap in units of i?Lei 
and Ca{9p, ijjs) is the arc length for the peak emissivity spot 
P{9p, i/jg) on this open field line. In principle, the peak posi- 
tion P{Op, ips) is dependent on the acceleration electric field 
and the emission mechanism. Based on our 1-D calculation 
of the acceleration field (see Figures |2) and later the emis- 
sivity (see Figure |8] later), the peak emission comes near the 
null charge surface. The height rp.AG for emission peak on 
open field lines can be related to the height of the null charge 
surface r-^{ips) by 



fp,AG(V's) = XKrj^ii>s) + (1 - A)KrN(O), 



(6) 



where k is a model parameter for the ratio of heights, and A 
is a model parameter describing the deformation of emission 
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Fig. 4. — The observed multi-band 7-ray light curves of the Vela pulsar (thick red lines in the right panels) and simulated profiles from our model (thin black 
lines). PI and P2 comes from the annular gap region, and P3 and bridge emission from the core gap region. The inclination angle a = 70° and the viewing 
angle ^ = 64° (Ng & Romani 2008) were used for modeling. The line of sight cuts across the 7-ray emission produced from only one pole. [See the electronic 
edition of the Journal for a color version of this figure.} 
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TABLE 1 

Model parameters for multi-band 7-ray light curves of the 
Vela pulsar 



GeV band 


K 


A 


e 


CTA 


^e, a" 


o'e. ci'' 


C2''' 


0.1-0.3 


0.68 


0.9 


1.17 


0.5 


0.0035 


0.0053 


0.009 


0.3-1.0 


0.70 


0.9 


1.20 


0.5 


0.0035 


0.0046 


0.009 


1.0-3.0 


0.72 


0.9 


1.15 


0.5 


0.0014 


0.0064 


0.01 


3.0-8.0 


0.72 


0.9 


0.88 


0.5 


0.0007 


0.0085 


0.006 


>8.0 


0.73 


0.9 


0.82 


0.1 


0.00085 


0.007 


0.003 



''The bunch scale for field lines in the annular gap. 
''The bunch scale for field lines of 90° < ips < 180° in the core gap. 
°The bunch scale for field fines of —180° < ipe < 90° in the core gap. 



location from a circle (see details in iLee et"an 120061 ). The 
emission peak position 'P' on each field line can be uniquely 

= ar csin [ ^ AG / , f ( a , -08 ) ] , where 



determined, i.e.. 



i?e,f (ct, Ips) is the field line constant of the open field line with 
■0s- Figure |3] shows the variations of rp AG (0s) and tn (V's) 
with t/'s. The minimum is at ips = 0° near the equator and the 
maximum at ■0s = ±180° near the rotation axis. 

The peak emissivity Ip{Op, ips) may follow another Gaus- 
sian distribution against d for a bunch of open field 
Unes (ICheng et al.ll2000l:lDyks & RudakllOOllFang & ZhangI 
l2OT0l) . i.e.. 



/p(6'p,0s) = /oexp 



OspiA) - 6*0? (■0s) 



(7) 



where /o is a scaled emissivity, ag is a bunch scale of 6 (in 
units of rad) for a set of field lines of the same 0s. 6'sp is 
used to label a field line in the pulsar annular regions, 6cp = 
(^n.Vb + ^p,iaJ/2 (i.e. rove (tps) = 0.5) is the central field 
line among those field lines with ip^. 

As seen above, we use two different Gaussian distributions 
to describe the emissivity on open field lines for both the an- 
nular gap and the core gap. The model parameters are inde- 
pendently adjusted to maximally fit the observed 7-ray light 
curves. In the core gap, we assume that the height of emission 
peak CG = e ■ ^p.AG, where e is a model parameter We 
adopted two different cre,c for the core gap because of the dif- 
ferent acceleration efficiencies for field lines in the two ranges 
of ^s. We will write ere. a for the annular region and ag^ c for 
the core region. 

3. To derive the "photon sky-map" in the observer 
frame, we first calculate the emission direction of each emis- 
sion spot ne in the magnetic frame; then use a transfor- 
mation matrix to transform into ngpin in the spin 
frame; finally use an aberration matrix to transform ngpin to 
= {hx, Hy, Hz} in the observer frame. Here 0o = 



^observi 

arctan(ny/nx) and ( ~ arccos(nz/ y^n^'^ + Uy^ + n^^) are 
the rotation phase (without retardation effect) of the emission 
spot with respect to the pulsar rotation axis and the viewing 
angle for a distant, nonrotating observer The det ailed cal- 
culatio ns for the aberration effect can be found in iLee et all 
(l2OT0h . 

4. We add the phase shift (50rct caused by the retarda- 
tion effect, so that the emission phase is = 0o — ^(pict- 
Here is no minus sign for 0o beacause of the different coor- 
dinate systems between our m odel and the outer gap model 
(iRomani & Yadigaroglull995h . 

5. The "photon sky-map", defined by the binned emission 
intensities on the (0, Q plane, can be plotted for 256 bins 
(see Figure lU. The corresponding light curves cut by a line 



of sight with a viewing angle ( = 64° are therefore finally 
obtained. For the viewing angle ( = 64°, any magnetic in- 
clination angles of a between 60° and 75° in the annular gap 
model can produce light curves with two sharp peaks and a 
large peak separation (e.g., 0.4 - 0.5), similar to the observed 
ones. The emission from the single pole is favored for the 
Vela pulsar in our model. 

The modeled light curves are presented in Figure 2] (black 
solid lines), with the model parameters listed in Table [T] 
Emission of PI and P2 comes from the annular gap region 
in the vicinity of the null charge surface, and P3 and bridge 
emission comes from the core gap region. The higher energy 
P3 emission (> 3 GeV) comes from lower height, whereas 
the lower energy 7-ray emission comes from a higher region. 
In the annular gap region, higher energy emission is mostly 
generated in higher region. Nevertheless, the 7-ray emission 
heights are above the lower bound o f the height determined 
by 7 - _B absorption (iLee et al.ll2010l) . 

The peak emission comes from different field lines and 
emission heights in the annular gap. The deformation of ra- 
diation beam i s related t o high value of geometric factor A 
as discussed in lDu et all (120101) . Owing to the aberration and 
retardation effects, the enhanced gamma-ray emission in the 
outer rim of photon sky-map make the peak very sharp, espe- 
cially for P2. For the Vela pulsar, the high inclination angle 
of about a = 70° is important to get the observed two sharp 
peaks with a large separation. 

3.2. Radio Lag 




Pulsar Phase 

Fig. 5. — The jointly modeled radio and 7-ray light curves (black solid 
lines) for the Vela pulsar. The radio emission is prodeuced from a higher and 
narrower region in the annular gap region of the same magnetic pole as the 7- 
ray. The observed radio profile (Ijlue dashed fine) for the Vela pulsar is taken 
from the website http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ephems/ and 
the 7-ray profile is observed by the Fermir[See the electronic edition of the 
Journal for a color version of this figure.} 
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With well-coord inated efforts for pulsar timing program, 
lAbdo et all ( 120 lOd) determined the phase lag between radio 
emission and 7-ray light curves. The radio pulse comes ear- 
lier by a phase of ^ 0.13 (see Figure|5]l. 

Radio emission might be generated in the two locations 
of a pulsar magnetosphere. One is the traditional low- 
height polar cap region for lo ng-period (P ^ Is) pulsars 
dRuderman & Sutherlandlll975l) . The other is the outer mag- 
netospheric region with high altitudes near the light cylinder 
dManche steij 120051) . For the polar cap region, the low radio 
emission height leads to a small beam, wh ich probably does 
not point to an observer for the Vela pulsar. iRavi et all (120101) 
propose that the radio emission from young pulsars is radi- 
ated in a high region close to the null-charge surface, i.e. the 
similar region for 7-ray emission. This is somehow similar 
to our annular gap model, in which the radio emission origi- 
nates from a higher and narrower region than that of the 7-ray 
emission. 

The modeled radio and 7-ray light curves in the two-pole 
annular gap model are shown in Figure |5] The region for 
the radio emission is mainly located at a height of ~ i?LC 
on certain filed lines with t/jg = —138°. Our scenario of ra- 
dio emission for the Vela pulsar is consiste nt with the narrow 
stream of hollow-cone-like radio emission (iDvks et al.ll2010l) . 
According to our model, not all 7-ray pulsars can be detected 
in the radio band, and not all radio pulsars can have a 7-ray 
beam towards us. 

3.3. 7-ray Spectra for the Vela Pulsar 

lAbdo et al] ( 1201 Obi) got high quality phase-resolved spec- 
tra (PI, P2, low-energy P3 and high-energy P3) and the 
phase-averaged spectrum of the Vela pulsar The observed 
7-ray emission is believed to originate from the curvature ra- 
diation of primary parti cles (iTang et al.l 120081: iHarding et alj 
l2008t IMeng eTani2008h . Here we us e the synchro-cu r vature 
radiation from primary particles (IZhang & ChengI 119951 
ICheng&Zhanglll996tlMengetal.ll2008l) and also the syn- 
chrotron radiation from secondary particles to calculate the 
7-ray phase-averaged and phase-resolved spectra of the Vela 
pulsar. 

We divide the annular gap region into 40 rings and 360 
equal intervals in the magnetic azimuth, i.e. in total 40x360 
small magnetic tubes. A small magnetic tube has a small 
area ^0 on the neutron star surface. From equation (|2}, the 
number density of primary particles at a height r is n(r) = 

^I'^jf^ cos (out, where c is the speed of light, and e is the elec- 
tric charge. The cross-section area of the magnetic tube at r is 
A{r) = BqAq/ B{r). Therefore, the flowing particle number 
at r in the magnetic tube is 



AiV(r)==A(r)As— ^cosCout, 
Zirce 



(8) 



here As is the arc length along the field. 

The accelerated particles are assumed to flow along a field 
line in a quasi-steady state. Using the calculated acceleration 
electric field shown in Figure |2] we can obtain the Lorentz 
factor 7 of the primary particle from the curvature radiation 
reaction 



7 = ( 



2e 



)3 =2.36x 1GV-5£;P-25 



(9) 



where pj is the curvature radius in units of 10^ cm and -Ey.e 
is the acceleration electric field in units of lO^Vcm"^. The 



pitch angle /3 of the prima ry particles flowing along a mag- 
netic field Une is dMeng et al. 2008) 



sin /? sa /3 w r/ 



7™cC 
e_B(r)p' 



(10) 



where 77 < 1, irie is the electron mass, and p is the curva- 
ture radius. The ch aracteristic energy jggy'i-cur synchro- 
curvature radiation (IZhang & Chenglll995l IMeng et al1l2008l) 
is given by 



2 p 



(11) 



X J(— + 1 -3-^)cos4/3 + 3-^cos2/3 + 4sin^/3: 

where tb = ^"'e'B(r)"^ '■^^ cyclotron radius of an electron, 
and h is the reduced Planck constant. 

The energy spectrum dA^/ d7 of t he a ccelerated primary 
particles is unknown. Harding et al.l (l2008i) have assumed it 
to follow a broken power-law distribution for pairs with in- 
dexes of —2.0 and —2.8 [see their equation (47)]. Here we 
assume the primary particles in the magnetic tube to follow 
one power law diV/d7 = Nqj^ with an index of F = —2.4. 
Here, A^o can be derived by integration the equation above us- 
ing the equations (HJ and (|9]l. The 7-ray s pectrum ernitted b y 
the primary particle can be calculated by (iMeng et al.ll2008l) 



F{E. 



dN 

d7 ^ 



-^E., (12) 



rc 




xG{x) 



■^^2/3(2;) 



d7, 



where Ail is the solid angle of the 7-ray beam, h is the Planck 



constant, x 



syn— cur 



+00 



^5/3 



[z)dz, 



K^/j,{z) and ^^2/3(0;) are the modified bessel function with 
the order of 5/3 and 2/3, and rc and Q2 are given by 



rc 



no 



[irB + p)nl- 



CCOSi 



eB{r) 
jrUcC ' 



TbX P'^ P TB j 

respectively. 

The secondary particles can be generated with a large mul- 
tiplicity (10'^ — 10"*) via the 7 — i? process in the lower regions 
of the annular gap and the core gap near the neutron star sur- 
face. Here, we assume that the energy spectrum of secondary 
particles follow a power-law, with an index of Fgcc = —2.8 
and a multiplicity of A/goc ~ 1000. The pitch angle of pairs 
increase due to t he cyclotron resonan t absorption of the low- 
energy photons (IHarding et al.ll20()8h . The mean pitch angle 
of secondary particles is about 0.06, adopted from equation 
( [Tol l with a slightly large factor ?/ > 1 owing to the effect 
of cyclotron resonant absorption. The synchrotron radiation 
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TABLE 2 

Best fit parameters for modeling the phase- averaged spectrum of the Vela pulsar. 







Rc 


r 


pri 
'mm 


pri 
'7m ax 


An 


2nd 
'mm 


2nd 

'max 


PI 


-no 


1.095 


0.62 


0.50x10^ 


1.79x10'^ 


0.11 


8.45x10^ 


1.31x106 


P2 


131 


1.278 


0.75 


0.30x10^ 


2.65x10'^ 


1.05 


6.25x10^ 


1.75x10*5 


P3 


-104 


1.122 


0.28 


0.35x10^ 


2.40x10'^ 


1.09 


3.05x10^ 


1.22x10^ 



Note. — Rc is the field line constant in units of -Rlc. and ^ is the emission height of a profile component. 




0.1 1 10 0.1 1 10 

Energy (GeV) Energy (GeV) 




Energy (GeV) Energy (GeV) 

Fig. 6. — Fitting the 7-ray phase-resolved spectra for the Vela pulsar. The observed data were taken from lAbdo et al.l fcOlObI) and plotted as black points with 
error-bars. The contributions of 7-ray emission from the both primary particles and secondary particles are plotted separately for profile components, PI, P2, 
P3 low-energy and P3 high-energy. The emission from the secondary particles only partly contribute the lower energy band. [See the electronic edition of the 
Journal for a color version of this figure.} 

ergy of < 50 GeV always have a t^-b <C 1 if the emission 
height is greater than a few hundred kilometers. 

To reduce the computation time, we calculate the synchro- 
curvature radiation at the "averaged emission-height" for 
three components, PI, P2 and P3, of the 7-ray light curve of 
the Vela pulsar. For PI, the emission height is about 0.62i?LC 
on the field line of a magnetic azimuth = —110°; for P2, 
the emission height is 0.75i?LC on the field line of if) = 131°; 
and for P3, the emission height is about 0.28i?LC on the field 



from secondary particles have some contributions to the low- 
energy 7-ray emission, e.g. < 0.3 GeV. 

We fur ther checked the optical depth t^-b of the 7— B ab- 
sorption (iLee et al.ll20Tol) 

1.55xl0V o ,2.76 X 10V/2pi/2 

here is in units of MeV, Bq, 12 is in units of 10^^ G. We 
found that the Fermi 7-photons of the Vela pulsar with an en- 
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Duet al. 2011 



T ' ' ' — I ' ' 

« Vela Fermi data 
_ . _ Calculated PI spectrum 




0.1 1 10 



Energy (GeV) 

Fig. 7. — Fitting the 7-ray pliase-averaged .spectram for the Vela pulsar. 
The observed data were taken from|Sbdo et al. (2010b) and plotted as black 
points with en'or-bars. The spectra for three profile components and the total 
phase-averaged spectram of the Vela pulsar are modeled from the synchro- 
curvature radiation from primary particles and synchrotron radiation from 
secondary particles. [See the electronic edition of the Journal for a color 
version of this figure.} 

line of t/j = —104°. We compute E\\ for the three peaks, and 
adjust the minimum and maximum Lorentz factor for primary 
particles, 7^^^'^^ and 7^,"^,' the minimum and maximum 
Lorentz factor for secondary particles, 7^^^"^^ and and 
the 7-ray beam angle Afi to fit the 7-ray spectra for the Vela 
pulsar 

We fitted the phase-averaged spectrum and phase-resolved 
(PI, P2, low-energy P3 and high-energy P3) spectra of the 
Vela pulsar as shown in Figure |6] and Figure |7] The best 
fit parameters for phase-resolved and phase-averaged spectra 
are similar as expected, and are listed in Table |2] The max- 
imum Lorentz factor of primary particles 7,p"x is consistent 
with that obtained from the curvature radia tion b alance of the 
outer magnetosphere models given by lAbdo et al. (2010bi)- 
The modeled spectra are not sensitive to 7^^"^^ or but 
quite sensitive to 7,p"x which is chosen around the value of 
the steady Lorentz factor given by equation (|9}. The solid 
angle AO was always assumed to be 1 by many authors for 
simplicity. We adjusted it as a free parameter around 1 for 
different phases. 

The synchro-curvature radiation from primary particles is 
the main origin of the observed 7-ray emission, while the syn- 
chrotron radiation from secondary particles can contribute to 
the lower energy band to improve the fitting. The peak ra- 
tios of PI and P2 shown in Figure |7] are roughly consistent 
with observations except for the band of 0.3-1 .0 GeV (cf. Fig- 
ure |4]i. The high-energy P3 is generated in the relatively low 
height of the core gap, where the particles have a higher ac- 
celeration efficiency than those for the low-energy P3, which 
lead to their cutoff energy different. The phase-resolved spec- 
tra for both high-energy P3 and low-energy P3 can be ex- 
plained in the synchro-curvature radiation from primary par- 
ticles from the core gap, with little contribution from the syn- 
chrotron radiation of secondary particles because they in gen- 
eral have small pitch angles with respect to field lines and 
large curvature radius. However, the synchrotron radiation 
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Fig. 8. — The emission at 1 3 and 8 GeV for the PI {top pannel) and 
P2 (bottom pannel) is not uniform along open field lines. It varies with the 
height. [See the electronic edition of tlw Journal for a color version of tliis 
figure.] 

from secondary particles does contribute to the 0.1 — 0.3 GeV 
band for PI and P2. 

In Figure [8] we plot the emission fluxes of PI and P2 com- 
ponent at 1, 3 and 8 GeV of the Vela pulsar against the emis- 
sion height. It is not uniform along an open field line. The 
bump at a low height for high energy 7-ray (e.g. >8 GeV) 
due mainly to the small curvature radius and large acceler- 
ation electric field there. In Section 3.1, we roughly took a 
Gaussian distribution along the arc (equation|5j to decribe the 
emissivity near the peak emission region, which is natural in 
our annular gap model and independent of the model paran- 
meters. 

4. DISCUSSIONS AND CONCLUSIONS 

The detailed features of 7-ray pulsed emission of the Vela 
pulsar observed by Fermi provide challenge to current emis- 
sion models for pulsars. 

The charged particles can not co-rotate with a neutron star 
near the light cylinder, and must flow out from the magne- 
tosphere. To keep the whole system charge-free, the neutron 
star surface must have the charged particles flowing into the 
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magnetosphere. We found that the acceleration electric field 
£■11 in a pulsar magnetosphere is strongly correlated with the 
GJ density pgj near the light cylinder radius Ri^c, while pgj 
at (i?Lc) is proportional to the local magnetic field Blc- It 
has been found that the Fermi 7-ray pulsars can be young pul- 
sars and millisecond pulsars which have a high Blc- This 
means that the acceleration electric field in a pulsar mag- 
netosphere is related to the observed Fermi 7-ray emission 
from pulsars. 

To well understand the multi-band pulsed 7-ray emission 
from pulsars, we considered the magnetic field configuration 
and 3-D global accleration electric field with proper boundary 
conditions for the annular gap and the core gap. We devel- 
oped the 3D annular gap model combined with a core gap to 
fit the 7-ray light curves and spectra. Our results can repro- 
duce the main observed features for the Vela pulsar. The emis- 
sion peaks, PI and P2, originate from the annular gap region, 
and the P3 and bridge emission comes from the core gap re- 
gion. The location and intensity of P3 are related to the emis- 
sion height in the core region. The higher energy emission 
(> 3GeV) comes from lower regions below the null charge 
surface, while the emission of lower energy of less than 3 GeV 
comes from the region near or above the null charge surface. 
Radio emission originates from a region, higher and narrower 
than those for the 7-ray emission, which explains the phase 
lag of ~ 0.13 prior to PI, consistent with the model proposed 



bv lDvks et aLrJIoTol) . 

Synchro-curvature radiation is a effective mechanism for 
charged particles to radiate in the generally curved magnetic 



field lines in pulsar m agnetosphere (IZhang & Cheng! ri995t 
ICheng & Zha"na[T996l) . The GeV band emission from pulsars 
is originated mainly from curvature radiation from primary 
particles, while synchrotron radiation from secondary parti- 
cles have some contributions to the low-energy 7-ray band 
(e.g., 0.1 — 0.3 GeV). Moreover, contributions of curvature 
radiation from secondary particles and inverse Compton scat- 
tering from both primary particles and secondary particles 
could be ignored in the 7-ray band. The synchro-curvature ra- 
diation from the primary particles and synchrotron radiation 
from secondary particles are calculated to model the phase- 
resolved spectra for PI, P2 and P3 of low-energy band and 
high-energy band and the total phase-averaged 7-ray spec- 
trum. 

In short, the 7-ray emission from the Vela pulsar can be 
well modeled with the annular gap and core gap. 
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